Deuteron production and elliptic flow in relativistic heavy ion collisions 



Yongseok Oh, x ' 2 'B Zi-Wci Lm. 3 ^ and Che Ming Ko 1 '! 

Cyclotron Institute and Physics Department, Texas A&M University, College Station, Texas 77843, USA 
Korea Institute of Science and Technology Information, Daejeon 305-806, Korea 
3 Department of Physics, East Carolina University, Greenville, North Carolina 27858, USA 

(Dated: November 19, 2009) 

The hadronic transport model ART is extended to include the production and annihilation of 
deuterons via the reactions BB <-> dM, where B and M stand for baryons and mesons, respectively, 
as well as their elastic scattering with mesons and baryons in the hadronic matter. This new 
hadronic transport model is then used to study the transverse momentum spectrum and elliptic 
flow of deuterons in relativistic heavy ion collisions, with the initial hadron distributions after 
hadronization of produced quark-gluon plasma taken from a blast wave model. The results are 
compared with those measured by the PHENIX and STAR Collaborations for Au+Au collisions at 
^/s NN = 200 GeV, and also with those obtained from the coalescence model based on freeze-out 
nucleons in the transport model. 

PACS numbers: 25.75.-q, 25.75.Ld, 25.75.Dw 



I. INTRODUCTION 

An important observable in heavy ion collisions at the 
Relativistic Heavy Ion Collider (RHIC) is the azimuthal 
anisotropy of the momentum distributions of produced 
particles in the plane perpendicular to the beam direc- 
tion, particularly the so-called elliptic flow (v 2 ) that cor- 
responds to the second Fourier coefficient in their az- 
imuthal angle distribution [J 0|. The measured ellip- 
tic flow is not only large but also shows a constituent 
quark number scaling, especially at intermediate trans- 
verse momenta, i.e., the dependence of the elliptic flows 
of identified hadrons on their transverse momenta be- 
comes similar if both are divided by the number of con- 
stituent quarks in a hadron. This scaling behavior of 
hadron elliptic flows is well described by the quark co- 
alescence model for hadron production from the quark- 
gluon plasma (QGP)formed in relativistic heavy-ion col- 
lisions @, 0, H, @, 0, Hi • Also, the measured elliptic flows 
of identified hadrons follow a mass ordering at low trans- 
verse momenta, namely, the strength of the elliptic flow 
becomes smaller as the hadron mass increases. This has 
also been well described by the transport model Q as 
well as by the ideal hydrodynamics [iJJ, . 

Recently, the elliptic flow of deuterons has been 
measured in Au+Au collisions at energy y/s NN = 
200 GeV [H [13 . The data from the PHENIX Col- 
laboration covers the intermediate transverse momentum 
(p T ) region [l2j |. while the STAR Collaboration has made 
measurements in a wider range of p T including the low 
p T region [l3|. The two measurements agree well at 
intermediate p T region (p T > 1.5 GeV/c) except that 
the STAR data show a negative elliptic flow at low p T 
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(p T < 1 GeV/c). A negative elliptic flow has also been 
seen in the preliminary data from the PHENIX Collab- 
oration for J/i/fs at p T ~ 1.5 GeV/c [2 EH- Negative 
values of deuteron elliptic flow cannot be explained by 
the coalescence model unless nucleons have negative el- 
liptic flows [HI , while negative J / if) flows may [l3] or may 
not [l8| require negative charm quark elliptic flows. Since 
measured nucleon elliptic flow does not show anynega- 
tive value in small transverse momentum region [12L 1 1 Qf | , 
observed negative values of deuteron elliptic flow raise an 
interesting question on the mechanism for their produc- 
tion and interactions in relativistic heavy ion collisions. 

In Ref. PH, two of us have studied deuteron produc- 
tion at RHIC in a dynamical model which is based on 
the time-dependent perturbation theory. Using the el- 
ementary reactions of NN dir, NNN — ► dN, and 
NNt: — > dir, we have computed the production rate of 
deuterons using the measured nucleon transverse momen- 
tum distribution that is parameterized by an effective 
temperature and a momentum-dependent elliptic flow. 
Although the energy is conserved in this approach, in 
contrast with the coalescence model, it needs the intro- 
duction of the reaction time and volume as parameters to 
fix the multiplicity of deuterons. The resulting deuteron 
p T spectrum and elliptic flow are found to be similar to 
those of the coalescence model based on the same nucleon 
momentum distribution. Compared to the experimental 
data [II, EH, the calculated deuteron spectrum is, 
however, too soft. As discussed in Ref. [la], this is due 
to the use of an effective temperature to model the effect 
of radial flow, resulting in the absence of correlations be- 
tween nucleon positions and momenta. For the elliptic 
flow of deuterons, this approach describes reasonably the 
experimental data except at small p T , where it gives pos- 
itive values while the data show negative values [l3j . 

To overcome the shortcomings of the model in Ref. [lj| > 
we use in the present work a transport model to study 
deuteron production and elliptic flow in relativistic heavy 
ion collisions. We assume that the particles produced 
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from hadronization of the created quark-gluon plasma 
are in thermal and chemical equilibrium and undergo col- 
lective motions, similar to those described by the blast 
wave model [20j , and they then interact with each other 
via hadronic rescattering. For the latter, we use a rela- 
tivists transport model (art) [2J [H, HH code embed- 
ded in a multiphase transport model (ampt) (24[. We 
also carry out a coalescence model calculation using the 
freeze-out nucleons to form dcuterons. This approach is 
different from the coalescence model used in Ref. [l6| as 
it includes the position and momentum correlations of fi- 
nal nucleons through their collective motions introduced 
in the initial blast wave model and further generated by 
final hadronic rescattering. 

This paper is organized as follows. In the next section, 
we describe the blast wave model used for generating the 
hadrons produced at hadronization that are taken as the 
initial conditions in the transport model for describing 
hadronic rescattering. Section HTII gives the results for the 
transverse momentum distributions and elliptic flows of 
pions and protons after hadronic scattering. The results 
for the deuteron transverse momentum distribution and 
elliptic flow arc presented and discussed in Sec. lIVl These 
results arc further compared in ScctionlVlwith those from 
the dynamical coalescence model. We also compare in 
Section [VTl the deuteron emission time distributions from 
these two models and discuss in Section IVlII the nucleon 
number of scaling of deuteron elliptic flow. Section I Villi 
contains a summary and discussions. The elementary 
production and annihilation processes for deuterons that 
are included in the transport model are described in the 
Appendix. 

II. INITIAL DISTRIBUTIONS OF HADRONS 

In relativistic heavy ion collisions, where a quark-gluon 
plasma is produced, hadrons are initially formed dur- 
ing hadronization of the QGP. To describe the distribu- 
tion of these hadrons, we use the blast-wave model which 
assumes that they are in thermal and chemical equilib- 
rium and undergo a collective expansion. As in the ART 
model, we include in the present work mesons such as 
7T, p, lo, r), K, K* , and 4> an d baryons such as N, A, 
A, and E as well as their anti-particles. We also include 
both deuterons and anti-deuterons in the initial distribu- 
tions. The Lorentz-invariant thermal distribution 
of these particles at temperature T c is then given by 

f(x,p)ocexp{-p> 1 u^/T c }, (1) 

where the four-momentum p 11 is 

P M = ijp iP) = ( m T cosh y,p T cos (f> p ,p T sin <f> p , m T sinhy) 

(2) 

and the flow four- velocity u^{x) is 

vf 1 = cosh ^(cosh ?7,tanhp cos </>, tanhp sin 0, sinhr/) 
= 7 r (cosh77,/3,sinh?7). (3) 



In the above, y is the energy-momentum rapidity, m T 
(= yj m 2 + pip) is the transverse mass with m being the 
mass of the considered particle, 7 T = — /3 2 , and 

77 (p) is the longitudinal (transverse) flow rapidity. We 
have followed the usual convention to define the z-axis 
along the beam direction and the x-axis in the reaction 
plane of the collision. The angles <f> p and <fi in the above 
are the azimuthal angles of the momentum and position 
vectors of a particle with respect to the a;-axis. For the 
longitudinal flow, it is assumed to be boost invariant so 
that the longitudinal flow rapidity is equal to the energy- 
momentum rapidity, i.e., rj = y. Wc further assume that 
the distribution in the energy-momentum rapidity is uni- 
form in the midrapidity. 

Because of non- vanishing elliptic flow v 2 in non-central 
heavy ion collisions, which is defined as 

/ Pl-P 2 y \ ^ 

2 \pI+pII 

with p x and p y being, respectively, the projections of the 
particle transverse momentum along the x and y axes 
in the transverse plane, the transverse flow velocity is 
anisotropic with respect to the azimuthal angle tp. Wc 
thus parameterize the transverse flow velocity as 

= 0(r)[l + e(p T )coa(2<p)]ri, (5) 

where n is the unit vector in the direction of their trans- 
verse flow velocity /3, which is taken to be normal to 
the surface of the hadronic system to be defined below. 
We have also introduced a p T -dependent coefficient e to 
model the saturation of the resulting v 2 at large p T as 
observed in experiments. Specifically, we parameterize e 
as 

e(p T ) = c t exp(-p T /c 2 ). (6) 
We further parameterize the radial flow velocity as 

where Rq is a parameter related to the transverse size 
of the initial hadron distribution in space and r is the 
distance of the particle from the origin of the coordinate 
system in the transverse plane. 

For the spatial distribution of the particles produced 
at hadronization, it is assumed to be inside a cylinder 
with its axis along the longitudinal direction and having 
an elliptic shape in the transverse plane, as only mini- 
mum biased collisions are considered in the present study. 
These particles are further assumed to be uniformly dis- 
tributed in the transverse plane. This means that the 
initial hadrons arc uniformly distributed in the spatial 
region of 
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In terms of the spatial elliptic anisotropy [9( 

the spatial region as given by Eq. ((5|) can be rewritten as 
r < R [1 + s 2 cos(2^)] , (10) 

leading to 

A = R (l + s 2 ), B = R (l-s 2 ). (11) 

For the coordinate system defined above, the initial shape 
of the hadronic system in non-central collisions then has 
s 2 < 0. 

We further introduce a formation tim e t for hadrons 
produced from hadronization as in Ref. [24j. Because of 
boost invariance, the time and the position in z-direction 
of an initial hadron are given by t = t cosh y and z = 
t sinhy, where y stands for rapidity. 

Since we do not consider the dynamics of partons in 
QGP before hadronization, the number of hadrons in the 
initial state cannot be determined a priori. We thus treat 
the initial number of positively charged pions N n + as a 
free parameter so that its number after hadronic rescat- 
tcring is fitted to the measured pion multiplicity. The 
number of other hadrons after hadronization is then de- 
termined from the assumption that they are in thermal 
and relative chemical equilibrium with respect to pions. 

For Au+Au collisions at center of mass energy 
y/s NN = 200 GcV, the parameters for the initial con- 
ditions are taken as follows. For the spatial distribution 
of initial hadrons, we use the parameters 

R = 5.0 fm, s 2 = -0.05, r = 2.5 fm/c, (12) 

similar to the quark distributions obtained in the AMPT 
model at hadronization. For the transverse momen- 
tum distribution, we use the critical temperature T c = 
175 MeV and 

/3 = 0.55, c x = 0.43, c 2 = 0.85 GeV/c (13) 

for the flow velocity parameterized in Eq. ([5]) . The num- 
ber of initial positively charged pions is taken to be 
= 100. This then leads to a proton number of 
N p = 2.5 in the thermal model with zero charge and 
baryon chemical potentials. These pions and protons as 
well as other hadrons are then uniformly distributed in 
the rapidity region \y\ < 2.0. The resulting initial trans- 
verse momentum spectra and elliptic flows of midrapidity 
(\y\ < 0.5) 7r + 's and protons, whose numbers are 71 and 
6.4, respectively after including those from decays of me- 
son and baryon resonances, are shown in Figs. [T]and[2] 
These results will be compared in the next section with 
those of the freeze-out pions and protons after hadronic 
scattering. 
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FIG. 1: (Color online) Transverse momentum spectra of pions 
and protons in midrapidity \y\ < 0.5. Data are from Ref. [26( |. 



III. PION AND NUCLEON TRANSVERSE 
MOMENTUM SPECTRA AND ELLIPTIC FLOWS 

We first discuss the transverse momentum spectra of 
pions and protons after hadronic evolution. The results 
in the midrapidity region \y\ < 0.5 are shown in Fig. [1] 
by the dashed and solid lines for pions and protons, re- 
spectively, and they are comparable to the experimental 
data measured by the PHENIX Collaboration [2tj . It is 
seen that the final pion and proton transverse momen- 
tum spectra are not very different from the initial ones 
shown by the dotted and dash-dotted lines, respectively. 
This can be understood from the fact that the effect of 
decreasing temperature of the hadronic matter due to 
expansion is compensated by an increase in the collec- 
tive flow velocity. A similar conclusion was drawn in 
Ref. [27j], where the effect of hadronic rescattering was 
studied through the hadronic transport model JAM [28|. 
We note that the number of pions and protons in midra- 
pidity after hadronic rescattering are N(ir + ) w 74 and 
N(p) « 5.4, respectively, including those from decays of 
baryon and meson resonances, and are only slightly dif- 
ferent from their initial values. 

In Fig. [21 the elliptic flow v 2 of pions and protons are 
shown by solid lines as functions of transverse momen- 
tum. It is seen that our model describes reasonably the 
experimental data from the PHENIX Collaboration ex- 
cept for p T > 2 GeV/c, where it overestimates the mea- 
sured pion elliptic flow. We have found that it is very 
difficult to find the parameters used for the elliptic flow 
in Eq. ((5]) that would describe both the pion and pro- 
ton v 2 in all p T region. This may not be surprising as 
the pion and proton elliptic flow at high p T follow the 
quark number scaling and are thus unlikely to be de- 



4 



0.2- 




0.1- 



p\ (GeV/c) 



p T (GeV/c) 



FIG. 2: (Color online) Elliptic flow v 2 of pions and protons 
for rapidities \y\ < 0.5. Data are from Ref. [in- 



scribed by the blast wave model. To overcome this may 
require the use of the blast wave model for the quark dis- 
tributions and the quark coalescence model to generate 
the initial meson and baryon distributions. In contrast to 
the transverse momentum spectra, the final pion and pro- 
ton elliptic flows are different from the initial ones shown 
by dashed lines. One sees that the pion elliptic flow is 
increased by scattering. For protons, their elliptic flow 
after scattering also increases at large p T but becomes 
smaller at low p T , reflecting an enhanced mass-ordering 
effect as a result of increasing radial flow. 



IV. DEUTERON TRANSVERSE MOMENTUM 
SPECTRUM AND ELLIPTIC FLOW 

In this section, we study the p T spectrum of deutcrons 
and their elliptic flow. As for other hadrons, we have 
assumed that they are initially produced at hadroniza- 
tion and then undergo hadronic scattering through which 
they can be annihilated and also reproduced. To include 
these effects, we have modified the art code embedded in 
the AMPT model as described in detail in the Appendix. 
The resulting deuteron p T spectrum is given in Fig.[3]by 
the solid line and is seen in relatively good agreement 
with the PHENIX data [HI]. The calculated deuteron 
number of about 0.046 in the midrapidity is, however, 
somewhat larger than the measured one. Compared with 
the deuteron initial p T spectrum shown by the dashed 
line, the final one has a larger inverse slope reflecting the 
further development of the radial flow of deuterons by 
hadronic rescattcring. 

Shown in Fig. [4] by the solid line is the calculated fi- 
nal deuteron elliptic flow, which is seen to describe rea- 
sonably the experimental data. The negative values of 
deuteron elliptic flow at low Pr (< 1 GeV/c) reported 
by the STAR Collaboration [13| are, however, not re- 
produced by the present model. We first note that the 
deuteron elliptic flow from the initial blast wave model 
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FIG. 3: (Color online) Transverse momentum spectrum of 
deuterons in midrapidity \y\ < 0.5. Dashed and solid lines 
are for initial and final deuterons, while the dash-dotted line 
is from the coalescence model. Data are from Ref. I29II. 



shown by the dashed line does have small negative val- 
ues at very low p T . This is due to following combined 
effects. Because of the larger flow velocity in the reac- 
tion plane than out of the reaction plane, deuterons of a 
given transverse momentum p T in the fire-cylinder frame 
have a smaller transverse momentum in the local frame 
if they move in the reaction plane than if they move out 
of the reaction plane. Since the number of deuterons 
with transverse momentum p' T in the local frame is pro- 
portional to p' T exp(— m' T /T), it is an increasing or de- 
creasing function of p' T for p' T < s/md,T or p' T > \fm^T. 
While the resulting deuteron v 2 is positive in the latter 
case as there are more deuterons of transverse momen- 
tum p T moving in the reaction plane than moving out of 
the reaction plane, it is negative in the former case. This 
effect is, however, reduced during subsequent hadronic 
scattering as a result of increasing radial flow velocity 
and additional production of deuterons. Indeed, the fi- 
nal deuteron v 2 at low p T is close to zero. Similarly, 
the deuteron elliptic flow at larger p T is reduced after 
hadronic evolution as shown in Fig. [5] 

In Fig. [5[ we show the numbers of protons and 
deuterons in the midrapidity \y\ < 0.5 as functions of 
time. Both are seen to decrease gradually with time and 
this is partially due to scattering into larger rapidity re- 
gions. 



V. THE DYNAMICAL COALESCENCE MODEL 

It is of interest to compare the results from the trans- 
port model with those from the coalescence model. In 
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FIG. 4: (Color online) Same as Fig. [3]for the deuteron elliptic 
flow. Data are from Refs. [l^. [ilj]. 
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FIG. 5: (Color online) Proton (scaled by a factor of 1/150) 
and deuteron numbers versus time in midrapidity \y\ < 0.5. 



the coalescence model, deuterons are produced by recom- 
bination of nucleons at freeze-out using the sudden ap- 
proximation. The momentum distribution of produced 
deuterons in this model is given by [3, l3(i| 



dpi 



dx\dp\ dx 2 dp 2 
f w {x^x' 2 ,p' l ,p' 2 )5^(p d -p x -p 2 llA,) 



with g = 3/4 from the consideration of the spins of nu- 
cleons and deuterons. In the above, d & N p /dx\dp\ and 
d 6 N n /dx 2 dp 2 are, respectively, the spatial and momen- 
tum distributions of protons and neutrons in the fire- 
cylinder frame and fw (x[ , x' 2 ; p[ ,p' 2 ) is the Wigncr func- 
tion of the deuteron in its rest frame. Assuming that the 
wave function of the deuteron is given by that of a har- 



monic oscillator, its Wigncr function is then 

fwix'^x'^p'^p'J = 8cxp(-p 2 cr 2 - x 2 /a 2 ) , (15) 
where a = 1/ ^JJSjJ with the reduced mass [i = m N /2 and 

x = x' 1 -x' 2 , P=\ (Pi ~ P2) > ( 16 ) 

with x\ and p\ being the positions and momenta of 
the coalescing proton and neutron in the center of mass 
frame of produced deuteron. The oscillator frequency 
is determined by the charge root-mean-square radius of 
the deuteron, yj (r 2 ) , = 1.96 fm [3(|, which leads to 
lj = 8.06 x 1CT 3 GeV. We note that the above deuteron 
Wigner function reproduces very well that obtained in 
Ref. using a more realistic deuteron wave function. 

The calculated results for the p T spectrum and elliptic 
flow of midrapidity (|y| < 0.5) deuterons from the coa- 
lescence model are shown by dash-dotted lines in Fig. [3] 
and Fig. 0] respectively. It is seen that except for small 
p T the deuteron p T spectrum in the coalescence model 
agrees almost perfectly with that in the transport model. 
The total deuteron number in this case is about 0.051 
and is only slightly larger than that from the transport 
model. For the deuteron elliptic flow v 2 , although the 
two models give similar values in the small p T region, 
the one from the coalescence model is larger at higher p T 
and the deviation between the two becomes larger as p T 
increases. As a result, the coalescence model does not 
describe the experimental data as well as the transport 
model. However, neither model could give a negative v 2 
for deuterons withp T < 1.0 GeV/c, contrary to that seen 
in the experimental data from the STAR Collaboration. 



VI. DEUTERON EMISSION TIME 
DISTRIBUTION 

We have examined the distribution of the deuteron 
emission times, which are the times for their last colli- 
sions in the transport model, and the results are shown 
in Fig. [5] The freeze-out times of nucleons are defined in 
the same way and their distribution is also given in Fig.[S] 
for comparison. In the coalescence model, a deuteron is 
formed from a pair of freeze-out or emitted proton and 
neutron, so its emission time is given by the latest time of 
the two emitted nucleons. Results given in Fig. [5] show 
that the distribution of deuteron emission times in the 
transport model is rather similar to that of nucleons. In 
the coalescence model, although the deuteron emission 
time distribution peaks at t ~ 10 fm/c, similar to that of 
nucleons, the two are otherwise quite different. In partic- 
ular, there are more deuterons emitted at later times in 
the transport model than in the coalescence model. We 
note that the non-smooth deuteron early emission time 
distribution in the transport model is due to the incom- 
plete destruction of initial deuterons, particularly those 
produced near the surface of the fire cylinder. Introduc- 
ing a diffused hadron distribution in the initial state may 
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FIG. 6: (Color online) Normalized freeze-out time distribu- 
tions of midrapidity protons (filled circles) and deuterons in 
the transport model (open circles) and in the coalescence 
model (open squares). 

help smooth the deuteron emission time distribution dur- 
ing earlier times. 

VII. NUCLEON NUMBER SCALING OF 
DEUTERON ELLIPTIC FLOW 
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FIG. 7: (Color online) Scaled elliptic flow in (a) transverse 
momentum p T and (b) transverse kinetic energy. The nucleon 
number n N is I for the proton and 2 for the deuteron. 

As in the case of the elliptic flows of identified hadrons, 
where a scaling according to the number of constituent 
quarks in a hadron has been observed [3l| . it is of in- 
terest to see if there is a similar scaling of the proton 
and deuteron elliptic flows according to their nucleon 
numbers. This is shown in Fig. [7] for the elliptic flow 
per nucleon as a function of momentum per nucleon 
[panel (a)] or as a function of transverse kinetic energy 
[.Ekt = (m 2 + Pt) 1 ^ 2 — Tn] per nucleon [panel (b)]. ft is 
seen that the nucleon number scaling of the proton and 



deuteron elliptic flows is not quite realized in the trans- 
port model, particularly at large p T . This scaling is, 
on the other hand, very well satisfied in the results from 
the coalescence model, although they do not describe the 
data as well as the transport model. 



VIII. SUMMARY AND DISCUSSIONS 

We have investigated deuteron production in heavy ion 
collisions based on a hadronic transport model. This 
is carried out by modifying the ART code in the AMPT 
model to include the production and annihilation of 
deuterons via the reactions BB <-> dM involving baryons 
and mesons as well as their clastic interactions with 
these hadrons. For the initial hadron distributions af- 
ter hadronization of the quark-gluon plasma, we use the 
thermal distribution based on a blast wave model. 

The initial hadrons are uniformly distributed in an el- 
liptic shape with a spatial elliptic anisotropy of —0.05. 
With the parameters such as the initial temperature as 
well as the transverse flow velocity and anisotropy fit- 
ted to reproduce the measured p T spectra and elliptic 
flows of protons, we have examined the production of 
deuterons in relativistic heavy ion collisions. The calcu- 
lated deuteron p T spectrum is comparable to, although 
somewhat softer than, the experimental data measured 
by the PHENIX Collaboration. The deuteron elliptic 
flow from the transport model also agrees reasonably 
with that measured by the PHENfX Collaboration. How- 
ever, the deuteron elliptic flow was found to be positive 
at small p T in contrary to the negative elliptic flow of 
deuterons observed by the STAR Collaboration. 

We have also compared our transport model results 
with those of the coalescence model. While the two give 
very similar deuteron p T spectra, the elliptic flows ob- 
tained in the two model are quite different, particularly 
at high p T . In addition, the coalescence model is found to 
give almost exact nucleon number scaling of the elliptic 
flow, while the transport model causes a deviation of the 
elliptic flow from the exact nucleon number scaling, which 
is, however, closer to the measured data. These results, 
including the non-negative values of the elliptic flow at 
very small p T , are not much changed if the cross sections 
for deuteron production and scattering are changed by a 
factor of two. 

Our study has raised two interesting questions on 
deuteron production in relativistic heavy ion collisions. 
One is the relation between the transport model and 
the coalescence model for deuteron production, and the 
other concerns the negative deuteron elliptic flow at 
low p T . Since the binding energy of deuteron is small 
(~ 2.2 MeV), the coalescence model has been considered 
an appropriate approach to describe deuteron production 
in nuclear reactions [32| . This is indeed the case in our 
study for the deuteron transverse momentum spectrum 
as the two models give very similar results. However, the 
deuteron elliptic flows from these two models are quite 
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different, with the transport model giving a better de- 
scription of the experimental data in comparison with 
the coalescence model. Part of this difference may be due 
to the neglect of rescattering of deuterons produced in 
the coalescence model. Although deuterons are produced 
from freeze-out nucleons in the coalescence model, they 
may undergo additional scattering in the hadronic matte, 
leading to a later emission time as results from the trans- 
port model have shown. As to the negative deuteron el- 
liptic flow observed in the preliminary experimental data, 
it is not seen in either the transport model or the coa- 
lescence model. The lack of negative deuteron elliptic 
flow in the coalescence model is expected as a negative 
deuteron elliptic flow requires a negative nucleon elliptic 
flow, which is not observed in the experimental data. In 
the transport model, a negative deuteron elliptic flow is 
possible if the final radial flow velocity is smaller than 
that reached in the present study, following the discus- 
sions in the end of Section lTVl Whether this can still lead 
to a good description of the measured deuteron trans- 
verse momentum spectrum needs to be checked. Further 
investigations both in theory and experiment are thus 
required to understand the elliptic flow of deuterons at 
low p T and to shed more light on the mechanisms for the 
production of low p T hadrons in relativistic heavy ion 
collisions. 
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APPENDIX: REACTIONS FOR DEUTERON 
PRODUCTION, ANNIHILATION, AND ELASTIC 
SCATTERING 

The ART code [2l|, [H, HH implemented in the AMPT 
code includes the interactions of 7r, K, r\, p, u>, <j), K* , N, 
A(1232), 5u(1535), Pn(1440) as well as their antiparti- 
cles. To extend the model to include deuteron production 
and annihilation, we modify the ART code to incorporate 
following reactions, 

BB' -> Aid, Md -> BB', (A.l) 

where M = ir, p, uj, 77, and B and B' stand for baryons 
N, A, Pn(1440), and Sii(1535). For the cross sections 
of the reactions BB' — > Md, we assume that their an- 
gular integrated mean squared matrix elements that are 
averaged over initial and summed over final spins and 
isospins are the same as that for the reaction TV N — > dir 
at same center of mass energies. The cross sections for 
the inverse reactions Md — > BB' are then determined 
from the detailed balance. Experimentally, the cross sec- 
tions for both the reaction pp — > dir + [3J, [34| and the 




Vs (GeV) 

FIG. 8: (Color online) Ex per imental data on the total cross 
sections for pp — > dir + [33l. |34| and for dir + — > pp [35l. |36|. |37| . 
The solid and dashed lines are the parametrization given in 
the text. 

reaction ir + d — ► pp [35l l36l. Wf\ have been measured, and 
the former can be parameterized as 

a{jyp^dir + ) = -^f{s\ (A.2) 

where p N and p n are, respectively, the magnitude of the 
three-momenta of initial and final particles in the center 
of mass frame. The function f{s), which is proportional 
to the angular integrated mean squared matrix elements 
that are summed over initial and final spins for the reac- 
tion pp — > ir + d, is given by 

f{s) = 26exp[-(s- 4.65) 2 /0.1] +4cxp[-(s-4.65)72] 
+ 0.28cxp[-(.s-6)7l0], (A.3) 

where the squared center of mass energy s is in units of 
GeV 2 and f{s) is in units of mb. For the inverse reac- 
tion dir + — > pp, its cross section is related to that for 
pp — > dir + via the detailed balance, i.e., a{dir + — ► pp) = 
(2p 2 v /3p 2 r ) cr(pp — » dir + ). These parameterizations are 
compared with the experimental data in Fig. [51 The cross 
sections for the isospin averaged reactions NN — ► dir and 
ird — ► NN can then be obtained from a(NN — > dw) = 
2 <j{pp — > dir + ) and a {dir — > NN) = a{dir + — > pp). 

In addition to the production and annihilation pro- 
cesses for deuterons, we also include their elastic scatter- 
ing with mesons M and baryons B, 

M + d^M + d, B + d^ B + d. (A.4) 

As for the deuteron production cross sections from 
baryon-baryon reactions, the cross sections for deuteron 
elastic scattering with baryons and mesons are obtained 
from those for deuteron elastic scattering with nucleons 
and pions, respectively, assuming that at same center of 
mass energies they have same angular integrated mean 
squared matrix elements that are averaged over initial 
and summed over final spins and isospins. 
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FIG. 9: (Color online) Experimental data from Refs. [38|, |39| . 
|40( | on the total elastic cross sections for Nd — > Nd. The solid 
line is the parametrization used in this work. 
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For the empirical nucleon-deuteron elastic scattering 
cross section [38|, HH, H(| , it can be parameterized by 



(T 



{Nd^Nd) = 2500cxp[-(s- 7.93) 2 /0.003] 
+ 300 exp " 



[-(s - 7.93) 2 /0.l] + 10, 
(A.5) 



where the cross section is in mb and s in GeV 2 . The 
comparison with available experimental data are given 
in Fig. 1 

The experimental data on the total cross sections for 
the nd elastic scattering are very scarce. Therefore, 
we have to rely on phenomcnological calculations. In 
Rcfs. [4l], [42|], based on a combined analysis of the re- 
actions pp — > pp, nd — > nd, and nd — > pp, the helicity 
amplitudes of these reactions are determined which then 
leads to the total elastic scattering cross sections for nd 
scattering. We use the results of Refs. [Il|, [42j, which can 
be parameterized as 



<j(nd^nd) = 63exp[-(s-4.67) 2 /0.15] 

+ 15 exp [-(s - 6.25) 2 /0.3] ,(A.6) 



FIG. 10: (Color online) Experimental data from Refs. [43l. 144 
I45l.l4r3.l47l] on the total elastic cross sections for nd — ► nd. The 
solid line is the parameterization used in this work. 



where s is in GeV 2 and the cross section is in mb. The 
comparisons with the results of Rcfs. [4l|, \^ and with 
the data are given in Fig. [TTJ] 
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